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An analytic study is presented of the E(£)e Jahn - Teller ( JT) polaron, 
consisting of a mobile e g electron linearly coupled to the local e g normal 
vibrations of a periodic array of octahedral complexes. Due to the linear 
coupling, the parity operator /C and the angular momentum operator J 
commute with the JT part and cause a twofold degeneracy of each JT 
eigenvalue. This degeneracy is lifted by the anisotropic hopping term. The 
Hamiltonian is then mapped onto a new Hilbert space, which is isomorphic 
to an eigenspace of J belonging to a fixed angular momentum eigenvalue 
j > 0. In this representation, the Hamiltonian depends explicitly on j and 
decomposes into a Holstein term and a residual JT interaction. While the 
ground state of the JT polaron is shown to belong to the sector j = 1/2, 
the Holstein polaron is obtained for the "unphysical" value j = 0. The 
new Hamiltonian is then subjected to a variational treatment, yielding the 
dispersion relations and effective masses for both kinds of polarons. The 
calculated polaron masses are in remarkably good agreement with recent 
quantum Monte Carlo data. The possible relevance of our results to the 
magnetoresistive manganite perovskites is briefly discussed. 

PACS numbers: 71.38.+i, 63.20.Pw, 72.80.Ga 
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I. INTRODUCTION 



The Jahn- Teller (JT) effect describes the interaction of lattice vibrational modes with 
orbitally degenerate electronic states and thus refers to a particular type of electron-phonon 
coupling. 1 Although this effect has proved indispensable for a proper understanding of the 
physics of a variety of systems, ranging from paramagnetic ions in nonmagnetic crystals 2 to 
structural phase transitions, 3 it is fair to say that its role in condensed-matter physics has 
been marginal for a long time. 

For nearly a decade, however, the significance of the JT effect is undergoing a profound 
change, triggered by the discovery of superconductivity in the fullerides 4 and of very large 
("colossal") magnetoresistance (CMR) in the manganite perovskites. 5 Because of their high 
symmetry, both classes of compounds fulfil the requirement for a JT interaction to occur, 
and numerous experiments seem to indicate that this is, in fact, the case. Manifestations 
of the JT effect in the fullerides have been reviewed by O'Brien and Chancey, 6 those in the 
manganites by Millis. 7 So far, however, there is no consensus as to the relative importance 
of the JT coupling in these materials. In their search of the origin of the CMR effect, 5 e.g., 
Millis et al. 8 argued that double exchange, 9 designed as a mechanism to induce ferromagnetic 
order in doped manganites, is not sufficient to account for the resistivity data and suggested 
that JT polaron formation is essential, whereas other authors 10 invoke ferromagnetic spin 
polarons to explain the effect. Problems of this kind could possibly be resolved by means 
of a detailed analytic theory of the JT polaron, yielding the (approximate) ground-state 
energy together with the corresponding eigenvector. Since the eigenstates of JT systems 
are vibronic in nature, 1 they may give rise to unexpected results for expectation values and 
correlation functions. 

As a first step in this direction, we study the JT polaron of symmetry type E ® e, which 
is most conveniently introduced by recalling some basic properties of Lai-^Ca^MnOs, a 
representative of the manganite family. 7 Each unit cell of the crystal contains an octahedral 
Mn0 6 complex and an average number of 4-x d electrons. Since the Hund's rule coupling 
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is believed to be very strong, the spins of all the d electrons are ferromagnetically aligned. 
Due to the crystal field produced by the oxygen ligands (point group Oh) the D state of the 
free Mn ion splits into a threefold degenerate t2 9 and a twofold degenerate e g level. Three 
of the electrons go into the tightly bound t2 9 orbitals forming a core spin of magnitude 3/2, 
while the remaining 1-x d electrons occupy the e g orbitals and are mobile. To study the 
formation of the polaron we start with pure CaMnC>3 (x — 1), where only the ti g orbitals 
are filled, and imagine that one additional electron is injected into the system (e.g., by 
replacing one Ca ion by La). The extra electron must go into the e g levels and, by virtue of 
symmetry, may couple to the e g normal vibrations of the octahedral complex. This type of 
vibronic interaction, where both the electron and the vibrational modes are of e g symmetry, 
is referred to as E <g> e JT coupling. In addition, the electron is allowed to move in a band 
composed of the local e g doublets. The resulting quasiparticle, consisting of the mobile e g 
electron and the concomitant e g distortion of the Mn0 6 octahedra, is designated as E ® e 
JT polaron. A somewhat simpler system, the E<S>b JT polaron, where b denotes a non- 
degenerate representation of the tetragonal site group, has already been treated by Hock, 
Nickisch and Thomas 11 nearly two decades ago. 

In Sec. II we introduce our model, together with the angular momentum operator J and 
the parity operator /C. Since only linear JT coupling is considered, these operators commute 
with the JT Hamiltonian, but not with the (anisotropic) hopping term. In Sec. Ill we 
show that K, generates new fermion operators such that both J and the JT term assume 
diagonal form with respect to the new fermionic basis. The spectrum of J is determined in 
Sec. IV, where we recover the well-known result that all eigenvalues of J are half-integral. 
We also show by rather general arguments that each eigenvalue of the JT Hamiltonian is 
still twofold degenerate. This degeneracy will be lifted by the hopping term. In Sec. V 
we construct a representation of the original Hamiltonian on a new Hilbert space, which is 
isomorphic to an eigenspace of J belonging to a fixed angular momentum eigenvalue j > 0. 
The new Hamiltonian depends explicitly on this quantum number and decomposes into a 
quasi-Holstein term and a residual JT interaction. This is the optimal form, which can be 
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reached by purely analytic means, and elucidates the close relationship between the JT and 
the Holstein polaron. While the latter is obtained for the "unphysical" value j — 0, the 
ground state of the JT polaron is shown to belong to the sector j = 1/2. A variational 
treatment of the new Hamiltonian is outlined in Sec. VI, where also some ground-state 
properties like the dispersion relations and effective masses for both kinds of polarons are 
presented. Our results are summarized in Sec. VII. 



II. THE MODEL HAMILTONIAN 

In the E (g) e JT polaron the state of the electron is completely specified by the vectors 
|ry), where i denotes the cell index and 7 = x, z the components of the e g doublet (spin 
indices are omitted since only a single e g electron is considered). The wave functions (f\ix) 
and (r\iz) transform like the orbitals d x 2_ y 2 and d 3z 2_ r 2, respectively, forming a local basis 
of the Eg representation associated with each unit cell. The state of the electron may be 
more conveniently specified by the operators ej 7 and e i7 , where ej 7 (e i7 ) creates (annihilates) 
an e g electron in the orbital state 7 = x or z at lattice site i. Similarly, the e g distortions 
of the MnOe octahedra may be described either in terms of the local normal coordinates 
Qi 7 or, more conveniently, by the bosonic creation and annihilation operators af 7 and a« 7 , 
respectively, where the indices have the same meaning as above. The Hamiltonian used for 
the description of the E ® e JT polaron reads 

TL = Tit + H v + Hjt, (1) 

where 

H t = -tJ2 ^ ■ h a ■ e t+a (2a) 

ia 

is the transfer or hopping term to be discussed below, 

n v = hnJ2^ - Si (2b) 

i 

describes the e g normal vibrations of frequency Q, while 
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Hjt = ghQ ^[( a lx + a ix)v x - (a\z + aiz)° z ]ei (2c) 

i 

represents the E ® e JT coupling, where a a (a = x, y, z) are the Pauli matrices and the 
coupling strength is expressed by the dimensionless parameter g. The column vectors 

e, = f 6%z j and a { = f j , (3) 

as well as their associated row vectors and aj, have been introduced for convenience and 
to avoid an accumulation of indices. 

The somewhat unusual form of the hopping term (2a), where the summation over a runs 
over the six nearest neighbors of site i, originates from the orbital degeneracy of the electronic 
states. In orbitally degenerate systems the transfer of electrons between neighboring sites 
depends on the orientation of the orbitals and the direction of the transfer. The matrix 
elements hi"'' of the matrices h a are entirely determined by symmetry. Their numerical 
values along the three cubic axes are tabulated in Ref. 12, where the interatomic matrix 
element Vdda is related to our hopping integral t. In the electronic basis defined by the vector 
e*j in Eqs. (3), the matrices h a take the explicit form 

\i ±x = {2a Q -V?>a x -a z )/A, 

h ±y = (2a° + V3a x -a z )/4, (4) 
h ±z = (a° + a z )/2, 

where a denotes the 2x2 unit matrix. 

In Eq. (2c) we have restricted ourselves to linear JT coupling, adopted by the majority 
of authors. 13 ' 14 Moreover, our model does not contain the intersite coupling of the normal 
modes, which one intuitively expects since oxygens are shared between adjacent MnC>6 
octahedra. These coupling terms will give rise to optical phonon branches and, as was 
recently pointed out by Hotta et al. 15 and Popovic and Satpathy, 16 to collective effects such 
as orbital ordering. Other terms like, e.g., the Hund's rule coupling should also be included 
in a more rigorous treatment. Hence, the neglect of all these couplings might render our 
model somewhat unrealistic; we hope to include these terms in future work. 
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Since Ti is invariant under lattice translations, the total crystal momentum P = K + Q 
is a conserved quantity. Here the operators K = X)k k & ' and Q = X) q Q d^ • <J q denote 
the crystal momenta of the electron and the e g vibrational modes, respectively, where e*k 
and a q are the Fourier transforms of e*j and <Jj. In addition to P, there are operators which 
commute with H v and Hjt, but not with the hopping term (unless TL t is assumed to be 
isotropic). Although these operators are not strictly conserved, they often greatly facilitate 
the diagonalization of H, as we shall see below. There are two operators of this kind which 
prove particularly useful: 
(i) the parity operator 



k. = gn, 



(5a) 



where 



and 



Q = exp ^7T ^ ' <*i 



1Z = exp 



(ii) the angular momentum operator 



where 



(5b) 



(5c) 



(6a) 



a J • o v ■ di 



(6b) 



is referred to as vibrational angular momentum. The spectral properties of these operators 
will be discussed in the following sections. Here it suffices to mention that quantities similar 
to /C and J also play an important role in isolated E(g)e JT centers, provided the JT 
coupling is linear. 17 In such systems the only eigenvalues of /C are 1 and —1, whereas those 
of J range over all half-odd integers. The most important properties of the operator /C may 
be summarized by the equations 
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lOeilC = a y ■ e*j, /C^fJj/C = —a*, (7) 

whose derivation rests on the well-known commutator expansion 

e s Ae~ s = A + [S, A] + (2!)- 1 [,S, [S, A]] + • • • . 

Using (7) we see that both Ti v + Hjt and J7" are left invariant by /C, i.e., W„ + Hjt, K- and 
J" form a complete set of commuting operators: 

[Hv + n JT , K] = [H v + Hjt, J] = [K, J] = 0. (8) 



III. GENERATION OF NEW FERMION OPERATORS 

In this section we shall exploit the properties of the parity operator /C to generate new 
fermionic creation and annihilation operators such that Hjt and J take diagonal form with 
respect to these operators. To this end we need another property of /C, which reads 

K? = 1 (9) 

and readily follows from Eqs. (7) and the relation (<r a ) 2 = a valid for all Pauli matrices. 
Hence, as in an isolated JT center, the only eigenvalues of K, are k — ±1. Moreover, since 
/C is also unitary, we have the additional relations /C = /C _1 = K) . 

As the next step, we need the projection operator V K for selecting the subspace associated 
with the eigenvalue k of /C. According to Lowdin, 18 V K is given by the expression 

V k = 1(1 + kK) (« = ±1), (10) 

which, apart from being Hermitian, has the properties 

E^ = 1 > ( lla ) 

K 

V K V K > = 5 KK ,V K . (lib) 

We also have the obvious relation JCV K = kV k , implying that the subspace projected out 
by V K is an eigenspace of K to the eigenvalue k. Property (lla) allows us to decompose 
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the Hamiltonian and the angular momentum operator J into components acting on the 
eigenspaces of /C as follows 

H = £ V K H t V K , +H V + J2 V K H JT V K , (12a) 

kk' k 

j = Y. v ^ v ^ ( 12b ) 

K 

where we have used Eq. (lib) and the fact that H v , Hjt, and J commute with V K . Since 
the hopping term does not commute with /C, the eigenspaces of the latter are mixed by Ht, 
as was to be expected. 

To obtain Eqs. (12) in explicit form, we need to calculate the operators e{P K . Using (7) 
and (10) we find 

^ * =5 (:;;«::)• < i3 > 

and we shall now prove that the products ICe^ (7 = x, z) on the right side of Eq. (13) may 
be replaced by Qe^, where Q is defined by Eq. (5b). To show this, let |\&) be an arbitrary 
vector of the underlying single-particle Hilbert space, 

|*> = £**r4r|0>> (14) 

where \I/j 7 are pure functions of the Bose operators a i7 , aj 7 and |0) denotes the common 
vacuum for all particles. If /Ce i7 is now applied to and use is made of the fact that 1Z of 
Eq. (5c) commutes with \I/j 7 (7 — x,z), the result is 

£e i7 |tt) = /C^ 7 |0) = GVyylO) = Ge^\m). 

Hence, JCe^ = Qe^ on the entire Hilbert space, which proves our claim. Expression (13) 
may thus be rewritten as 



eiP K = - = (15a) 

2 \e ix + iKQe iz 



where 



^fc) (15b) 



is a normalized vector (i.e., uj ■ u K = 1), while the quantities 

d iK = ~^( e iz ~ iKQe ix ) (15c) 

behave like ordinary fermion operators, i.e., 

[diK,dj K ,}+ = 5 i:j 5 KK ', [d iK ,d jK >] + = 0. (16) 

However, due to the presence of the operator Q in Eq. (15c), the di K , d} K cease to commute 
with the Bose operators a i7 and aj 7 , but continue to commute with quadratic forms like Q 
and the vibrational angular momentum M. of Eq. (6b). 

With the help of Eqs. (15), the various parts of the Hamiltonian (12a) may now be 
expressed in terms of the new fermion operators d iK and d} K {k = ±1). We start with the 
hopping term Tit, which is transformed into 

n t = -tJ2T,dLrr'd l+a , K/ , (17a) 

ia kk' 

rf =uJ-h a -u K ,, (17b) 

where the matrices h a are given by Eqs. (4). In the basis defined by the vector di = (* + ), 
Eqs. (17) may also be written as 

Ht = -t di ] ■ T a ■ di+a, (18a) 

ia 

where the new hopping matrices r a read 

T±x = (2a -o*- VZGa y )/4, 

T ±y = (2a° -a x + V3Qa y )/A, (18b) 
T ±z = (a° + <T x )/2, 

being now explicit functions of the operator Q. The vibrational part 7i v remains unchanged, 
while the JT coupling takes the form 

9 



Hjt = ghQ J2 d Li iK ( a lx + a ix)Q - (a\ z + a iz )]d iK . (19) 

in 

In deriving this result we have used the relations a i7 M K = -u_ K a i7 , uj ■ a x • U- K = —ikQ, and 
uj ■ a z ■ u_ K = 1. Finally, the angular momentum operator J of Eq. (12b) is obtained as 

J = M-{Q/2)Y, K did iK , (20a) 

in 

and we see that Tijx and J are now diagonal with respect to the new fermion operators. 
The vibrational part of J is, however, still nondiagonal. The diagonalization of M. is readily 
accomplished by means of the substitutions 

— ( 

&ix * (fliz &ix)i 

which leave H t and H v invariant, while M. is brought to the diagonal form 

M = Yl( a lx a ix ~ a\ z a iz ) = Mi. (20b) 

i i 

The JT coupling is transformed into the expression 

Hjt = -V2ghQ J2 dl(a tx U K + a iz Tl_ K + H.c.)d tK , (21) 

in 

where the new projection operators 

U k = ^(1 + kG) (« = ±1) (22) 

have been introduced. For later purposes we need the properties 

n K a^ = ai^H— K , (23a) 
ILJLy = 5 KK ,I1 K , (23b) 

U K + n_ K = 1, (23c) 

n K - n_ K = kg, (23d) 

which follow immediately from definition (22). Before we set out to develop strategies for 
dealing with the complicated vibrational terms in T~Cjt, we shall first derive the spectrum of 
J and investigate the possible symmetries and degeneracies of the JT Hamiltonian. 
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IV. SYMMETRIES AND DEGENERACIES 



In isolated JT centers there often exist "hidden" symmetries giving rise to unexpected 
degeneracies of the energy levels. In the linearly coupled E ® e JT center, e.g., all eigenvalues 
are twofold degenerate, and higher-order coupling terms are necessary to (partially) remove 
the degeneracy. 1 This section is devoted to a study of these degeneracies in JT crystals like 
the manganites. 

To investigate the possible symmetries of our system, we first need the spectrum of the 
angular momentum operator J , Eqs. (20). Since J , /C, and 7i v commute with each other, 
these operators possess common eigenstates. The latter are of the form 

K K ) = E^4I1K^) (« = ±1), (24a) 

i i 

where C iK are coefficients and the product extends over the local eigenstates of the isotropic 
oscillator in two dimensions: 19 

(J \(n i +m i )/2( n \ \{m-mi)/2 

Krm) = { lx \ 1 lz) |0), (24b) 



V v — 2 — ; ! 



rii = 0, 1, 2, • • • ; m i = n i ,n i -2,---, -rij. (24c) 

In fact, a simple calculation shows that 1^°^) is an eigenstate of /C, J, and H v with the 
respective eigenvalues k, j K , and E^ K , where 

E^ K = nhn, ^ = m -^(-l)» (25a) 



n 



= J2m, m = ^m i . (25b) 



There are many other linearly independent eigenvectors of /C, J , and H v belonging to the 
same eigenvalues: all vectors of the form (24a), whose quantum numbers rii and rrii satisfy 
the constraints (24c) and (25b), are also eigenstates with the required properties. Together 
they span a vector space U®, and we see that the energies E^ K are highly degenerate. 20 
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For a single site we have the relation rrii = rii — 2pi (pi = 0, 1, • • • , n») which, after 
summation over all cells, becomes m — n — 2p {p — 0, 1, • • • , n). Thus, for a given n, the 
quantum number m may take the n + 1 integral values m — n, n — 2, • • • , —n, whence we 
conclude that the angular momentum quantum numbers j K must all be half-integral, as in 
an isolated E ® e JT center with linear coupling. 1 A more detailed analysis of j K , Eq. (25a), 
requires a distinction between even and odd n (n and m always have the same parity, both 
being either even or odd). It is not difficult to verify that, for fixed k, both cases yield the 
same eigenvalues j K so that we may restrict ourselves to even m (n). The angular momentum 
quantum numbers may thus also be written as 

j K = m- k/2 (m = 0, ±2, ±4, • • •) 

or, explicitly: 

••.,-5/2,-1/2,3/2,7/2,... k = 1 
3k = < (26) 
•••,-7/2,-3/2,1/2,5/2,... k = 1. 

We now set out to examine more closely the structure of J and Hjt which, for the 

present purpose, are written as J = Y.kJ^ an d 'Hjt = Sk^jt- A glance at Eqs. (20) 

and (21) then reveals that the substitution d iK — > di- K , combined with the interchange 

o-ix Ojz (these operations correspond to canonical transformations and, hence, do not 

affect the eigenvalues), has the effect that Hj K T — > Hj~ T ^ and — > —J^ K \ Denoting the 

eigenvalues of Hj K T and by E n j K and j K , respectively, and the common eigenvectors of 

these operators by \^f n j K ), we may thus draw the following conclusions: 

1. Hjt an d H.jt have the same eigenvalues which must, therefore, be independent of k: 
E njK = E nj . Since the corresponding eigenvectors \^ n j+) and \^ n j-) are orthogonal by 
virtue of Eq. (lib), each eigenvalue E n j of Hjt is necessarily twofold degenerate. 

2. The spectra of and have the property that to any positive eigenvalue j + of 

there is always a negative eigenvalue j- = —j + of and vice versa. This property 
is most clearly reflected by Eq. (26). Since the eigenvalues of Hjt are independent of k, 
they can only depend on j = \j K \ = 1/2, 3/2, 5/2, • • •. 
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Thus, we see that the huge degeneracy of the eigenvalues E® jK of H v (see Ref. 20) is nearly 
completely lifted by the JT Hamiltonian. The remaining twofold degeneracy of H JT) which 
is of the same origin as that in the linearly coupled E ® e JT center, will be removed by the 
hopping term Tit (apart from accidental degeneracy). 

Only relatively few of the vectors contained in are simultaneous eigenstates of J 
and TCjt- To find these eigenvectors, we shall take advantage of the existence of a simple 
operator C, which also commutes with J and is of help to select the proper candidates. This 
operator will be shown in Sec. V to emerge from the JT Hamiltonian and reads 

C = £4(l-«0MKe. (27) 

in 

There are, in fact, two orthogonal sets of eigenstates both belonging to the same quantum 
number j, but to different eigenvalues of C. For positive j the two sets are represented by 
the vectors (similar vectors have been constructed in Ref. 21) 

oo 

=E E C+( ni )di\m + 2n i ,m)ni0i,0,>, (28a) 

i n;=0 tyi 
oo 

l*^> = E E Cf K {m)dl\m + 2n t + 1, m + 1) J[ |0,, t >, (28b) 

i rii=0 l^i 

where m is the same for all i and may assume the values 

m = j-l/2 = 0,l,2,.... (28c) 

It is then straightforward to verify the eigenvalue equations 

J\*%J=j\*ijJ, (29a) 
C\*% K } = (l/2±m% K ), (29b) 

whose validity requires that 

K = -(-l) m = (-\y +1 / 2 . (29c) 
Common eigenstates of J and C for negative j also exist, but are not needed here. 
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Hence, the vectors 1^^) are common eigenstates of J and C. Physically they represent 
polaronic states, where the electron is accompanied by an on-site distortion (vibrational 
excitation) of the molecular complex, all other complexes not coinciding with the location of 
the electron being left in their vibrational ground states. All vectors (28a), where m is given 
and k is fixed by Eq. (29c), form a subspace Uj~, while those obtained from span a 

subspace U~, which is orthogonal to Uj~ . The direct sum of these spaces will be denoted as 
Uj to remind us that this is an eigenspace of J to the eigenvalue j = m+ l/2>0. Since 
this is the most general eigenspace, which is compatible with the existence of the operator 
C, the simultaneous eigenvectors l^njn) of Hjt and J are necessarily all contained in Uj. 

V. REPRESENTATION OF U FOR FIXED QUANTUM NUMBER j 

Our main goal in this section is to construct a representation of the Hamiltonian on a 
subspace Vj, which is defined to be isomorphic to the space Uj introduced at the end of the 
preceding section. We start with the JT term, whose representation rests on the operators 

A= V2J2 d l(a ix Tl K + a iz Tl_ K )d iK , (30a) 

IK 

A ] = V2j2dL(a\ x U_ K + a\ z U K )d lK , (30b) 

IK 

allowing the JT Hamiltonian (21) to be written in the simple form 

H JT = -ghQ(A + A r ). (31) 

The merits of this seemingly trivial reformulation will become obvious later on. As the next 
step, we calculate the products AA^ and A^A. Using properties (23a) and (23b) we find 
the expressions 

AA ] = 2J2 <&(1 + oL««n_ K + a\ z a iz n K )d iK , (32a) 

in 

AU = 2J2 dL(a\ x a ix U K + a\ z a tz Il. K )d tK) (32b) 

in 

which enable us to set up the commutator [A, A' ] and the anticommutator [^4, A^ ]+. Using 
properties (23c) and (23d) we obtain 
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\[A,A^\ = £4(1 - KGMJdi* = C, (33a) 

in 

]-[A, ^ ]+ = E 4(1 + ^ ■ di)d lK = Af, (33b) 

where C = \[A, A^\ is the operator (27), whose eigenvalue problem is given by Eq. (29b). 
The result may be restated as follows: on the subspace Uj~ the operator C reduces to a 
positive integer and takes the form 

\[A,A^} =C = m+l=j + l/2 (m>0), (34) 

whereas on UJ it reduces to the number 1/2 — j, which is negative for all j > 1/2. To 
appreciate this result, we now investigate the anticommutator Af, Eq. (33b). First of all 
one realizes that Af is a positive operator. This property, together with the identity 

Af = A ] A + C, (35a) 

imposes a constraint on C requiring that this must also be a positive operator. The only way 
to guarantee that C and, hence, Af are positive for all j > 1/2 is to restrict the commutator 
to the subspace as in Eq. (34). 

Particularly interesting are the commutators [Af, A] and [Af, A^], since they agree with 
those of ordinary Bose operators. 22 In fact, using Eqs. (30) and (33b) one obtains 

[Af,A] = -A, [AT, A ] } = Al (35b) 

We also need to express TC V in terms of the operators A and A. To this end we calculate the 
expectation value of Af in the state \in) = d- K \0) e , where |0) e denotes the electronic vacuum. 
The result is (m\Af\iK) — 1 + • a j which, after summation over all sites, yields 

H v = HQ ]r((m|A/lm) - l). (36) 

i 

We are now in a position to construct the desired representation of the JT Hamiltonian, 
as defined at the beginning of this section. Our idea is to map the Hilbert space Uj onto a 
new space Vj, isomorphic to Uj, and to construct operators on the new space satisfying the 
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same relations as A and A^ on Uj (this procedure is closely related to that employed in the 
bosonization of spin operators). The new space Vj is spanned by all vectors of the form 

l*W = E E ^(^4-^=10), (37a) 
i r H =o V(2ni)! 

= E E c£fo)4 ; j = 10), (37b) 

i n i= ^(2^+1)! 

where b} and 4 create new bosons and fermions, respectively, and |0) is the common vacuum 
of the new particles. The new fermion operators are similarly defined as our previous d iK of 
Eq. (15c): 

c in = ( e * z — "i K G&ix)i (38a) 
G = exp (in b l b ^j ■ ( 38b ) 
The required isomorphism between Uj and Vj is achieved by the mapping prescription 

- (39) 

which obviously establishes a one-to-one correspondence between all vectors of Uj and those 
of Vj. As in our construction of the space Uj, the new space may also be conceived as the 
direct sum of two orthogonal subspaces, Vf and Vj , spanned by all vectors of the form 
(37a) and (37b), respectively. The quantum number n in Eqs. (37) may be assigned the 
value (29c), but we shall see below that the results are independent of this choice. 
On the new space Vj we now define the operator 

^. = £4A°'W, (40a) 

ire 

where is self-adjoint and explicitly dependent on the quantum number j: 

4^P- + (1 + 7 ^) 1/2 P + , (40b) 

P± = ±(1±G). (40c) 
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The operators P±, the analogues of our previous II K , satisfy a set of relations differing from 
Eqs. (23) only in notation (we have, e.g., that P±bi = 6jP T ). We now set out to prove that 
Aj and Aj satisfy the same algebraic relations on Vj as A and A^ on the space Uj (another 
proof of the equivalence of the operators A and Aj in terms of their matrix elements is given 
in Appendix A). To this end we start by calculating the products AjAj and AjAj. Using 
the well-known relation F(b}bi)bi = biF{b}bi — 1), we find that 

AjA) = £ 4(6/6, + 2jP + + l)c iK , (41a) 

in 

A)A 3 = £ clibjh + 2jP_)c iK . (41b) 

in 

By means of these relations, we may now evaluate the commutator [A,-, Aj] and the anti- 
commutator [Aj, Aj] + . A simple calculation gives 

\[Aj,A)\ = E 4(V2 + jG)c iK = Cj, (42a) 

%K 

\[Aj, A) ] + = £ 4(6/6, + j + l/2)c m ^ iV,, (42b) 

IK 

where Cj and Aj denote the analogues of our previous operators C and Af . 

If our formalism is to be meaningful, we expect the commutators [Nj, Aj] and [Aj, A-] 
to be the same as in Eqs. (35b), differing from the latter only in notation. This is in fact 
the case, for a straightforward calculation based on Eqs. (40a) and (42b) shows that 

[N j ,A j ] = -A j , [Nj,A)]=A], (43) 

and we recover the Bose-like commutation relations (35b). 22 It still remains to be verified 
that the commutator (42a) agrees with that of Eq. (34). This is readily shown and follows 
from the observation that on the operator G acts like the unit operator, while on Vj 
it has the eigenvalue —1. Moreover, since Nj is again a positive operator, Cj must also be 
positive by the same arguments as those used in conjunction with Eq. (34). Hence, on the 
allowed subspace Vf the commutator reduces to the positive integer 

l[Aj,A]]=Cj=j + l/2, (44) 
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exactly like C on Wj (see Eq. (34)). Thus, we have shown that the operators Aj, Nj and 
A, A\ Af are defined on isomorphic Hilbert spaces and satisfy the same algebra. The two 
sets are, therefore, physically indistinguishable from each other. In retrospect we realize that 
relations (43) and (44) do not depend on the sign of k in Eqs. (37). This independence is 
important, since it restores the twofold degeneracy of the eigenvalues of the JT Hamiltonian. 

The equivalence of the sets {Aj, Aj, Nj} and {A, A\ Af} will now be exploited to 
construct the desired representation of the Hamiltonian on the space Vj. Consider first the 
vibrational term TC V , whose representation on Uj is given by Eq. (36). To represent TC V on Vj, 
we start by using the correspondence Af ~ Nj and d iK ~ q k . Thereby the matrix element 
{in\Af\in) is mapped on (in\Nj\in), where \in) = cfjO) e and |0) e denotes the electronic 
vacuum in Vj. The representation of 7i v on Vj is then obtained as 

H®/(m) = EKwcMwe) - 1] = J2(blh+J - 1/2), (45) 

i i 

where Eq. (42b) has been used. To find the representation of the JT term, we must go back 
to Eq. (31) and exploit the correspondence A ~ Aj. In this way the representation of TijT 
on Vj is found to be given by the expression 

Hf T /{hV) = -g(A j + A]) = -sE4(4°\ + H.c.)q k , (46) 

in 

where use has been made of Eq. (40a). For later purposes it proves more convenient to 
rewrite Hj^ in the form 

Hf T = - ghn V cl(P + Bl% + R.c.)c iK , (47) 



where 



H% = -ghn X: 4(6/ + h)c iK (48a) 



and 



*M 1 + ;K7lf " L (48b) 
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We mention in passing that expressions quite similar to those in Eqs. (45), (47) and (48) 
have been derived for the isolated E ® e JT center by means of a rather different approach. 23 
In the unphysical limit j — the JT interaction reduces to the displaced-oscillator term H^, 
since B\ j) = in this case. The hopping term deserves some comment. 

In Sec. II it was pointed out that H t does not, in general, commute with J implying 
that Ht may induce transitions between subspaces belonging to different eigenvalues j of J . 
However, from the vibrational term (45) it follows that the ground state belongs to j = 1/2 
and is separated from the state with j = 3/2 by the energy NhQ, which tends to infinity 
in the thermodynamic limit. This gives rise to a kind of selection rule, allowing transitions 
only within the ground-state manifold, and we may, therefore, regard H t as being restricted 
to the subspace Uj=i/ 2 . Hence, its representation on Vj=i/2 takes the form 

H t = -tJ2& -ta-Z+a, (49) 

ia 

where the matrices t a differ from the T a in Eqs. (18b) only in the replacement of Q by G. 

Putting together our findings, the representation of the total Hamiltonian can now be 
easily written down. For the ground state (j = 1/2), the result may be recast into the form 

H = H QH - ghQ cl(P+B^bi + H.c.)c m , (50a) 

in 

where B\ j) (j = 1/2) is defined by Eq. (48b) and 

Hqh = -t XI Cj f • t • c i+a + httJ2 bjbi - ghVt J2 c\ K (b} + bi)c iK (50b) 

ia i in 

will be referred to as quasi-Holstein model. Although Hqh and the standard Holstein 
model 24 have the same formal appearance, they differ from each other in two respects. First 
of all they differ in the hopping term, which is isotropic in the standard model and anisotropic 
in our case. A more subtle difference lies in the properties of the fermion operators: while 
those of the standard model commute with the boson operators, this is not the case with 
the Qkj c\ k in Eqs. (50). The peculiar behavior of our fermion operators will, however, not 
entail any problems in the further analysis. 
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VI. VARIATIONAL TREATMENT 



The structure of the vibrational term (45) led us to conclude that the ground state of 
the E<S>e JT polaron belongs to j = 1/2 and thus agrees with the well-known fact that 
the lowest state of the isolated E(g)e JT center belongs to the same quantum number. 1 
Subsequently our Hamiltonian H of Eqs. (50) will be subjected to a variational treatment, 
yielding the (approximate) ground-state energy and the corresponding eigenvector for both 
the JT and the Holstein polaron. From Eqs. (50) we see that the quasi-Holstein model 
results for j = 0, while the JT case is obtained for j = 1/2, allowing both types of polarons 
to be treated in a unified way. A substantial simplification of the Hamiltonian may be 
achieved by an expansion of the square root in Eq. (48b) which, for j = 1/2, converges for 
all eigenvalues of b-bi. Keeping only terms linear in j, the expansion of becomes 



Since B\ 3 ' vanishes for j = 0, we expect Eq. (51) to be reasonable for j = 1/2. Indeed, for 
the isolated E ® e JT center this turns out to be an excellent approximation, 23 and we shall 
use it in the subsequent analysis. 

One of the main ingredients of our treatment is a unitary operator (Jost transformation) 
having the property that the transformed Hamiltonian assumes diagonal form with respect 
to the fermion momenta. This transformation, which has already proved its utility in various 
other polaron problems, 25 has the form 



where the last equality only holds on the single-particle Hilbert space. The operator 



Bf =j(b}b l + ^ + 0{f). 



(51) 




(52) 



Ui = exp(-iQ • RJ, 



(53a) 



where 



(53b) 



q 
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denotes the crystal momentum of the boson field, commutes with all q k and is a translation 
operator for the Bose particles: 

UfbtUi = &,_,. (53c) 

To better understand why the Jost transformation leads to a partial diagonalization of 
the Hamiltonian in momentum space, we recall from Sec. II that the total crystal momentum 
P = K + Q, where K = J2k k c\J ■ Ck denotes the fermion momentum and Q is given by Eq. 
(53b), is a conserved quantity. Using the alternative form U = Sk*^ • c*k+Q, one readily 
establishes the relation U'PU = K, showing that in the transformed system K plays the role 
of the total crystal momentum. This necessarily implies that the transformed Hamiltonian 
becomes diagonal in the fermion momenta k, the eigenvalues of K, whereas the indices k 
remain unaffected. 

To obtain the transformed Hamiltonian U^HU, we use the last equality in Eq. (52) and 
property (53c). The calculation is straightforward and leads to the result 

U ^HU = ~tJ2 U aCi ] ■ t a ■ c i+a + hilJ2 b}bi 

ia i 

-<7M2£4[(1 + P + Bi J % + H.c.]c m , (54) 

in 

where bo = 6j=o and is given by Eq. (51). After Fourier transformation of the fermion 
operators the Hamiltonian assumes the expected form 

^ff^EEiC'^, (55a) 

k kk' 

where 

= -T^ Q + 5 KK ,hn{J2b}b t 

i 

-g[(l + P + B^)b + R.c.]} (55b) 

and 

^k-q =«ECexp[i(k-Q) -RJ. (55c) 

a 
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Our variational treatment is based on the displacement transformation 

^ = E c Lj4c kK , (56a) 



.26 



defined on the single-particle Hilbert space, where the operator 



V k = exp 



N q 



(56b) 



displaces the 6 q according to the familiar formula V^bqV^ = 6 q + v k / x/iV (iV denotes the 
total number of sites). The unknown parameters v k will be determined from the variational 
principle. The Hamiltonian then acquires the form 

H = (UVyH(UV) = E E cL^k"V k tw, (57) 

k kk' 

where H£ K ' is given by Eq. (55b). Our further strategy may now be outlined as follows: as 
the first step, we introduce new fermion operators / kK , fl K by means of the relation 

/LlO) = (cV • L k ) K |0) = E^cLlO), (58) 

k' 

where L k is a unitary 2x2 matrix, which will be specified below. In terms of the new 
fermionic basis the Hamiltonian (57) takes the form 

H = E E /L(Lk^HkVkL k )^/w (59) 

k kk' 

and possesses the matrix elements 

(0|/ k ^/L'|0) = (4M k L k ) KK , = E^kT^'L^', (60a) 

w 

where the expectation values 

M£ X ' = (0|\/ k f ^ A V k |0) (A, A' = ±1) (60b) 
will be functionals of i> k . The matrix L k is now fixed by the requirement 

(4M k L k ) KK , = £ K (k)<W, (61) 
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which gives rise to two energy bands £ K (k) (k = ±1). The variational parameters i> k are 
then obtained from the condition 5£_(k)/5t> k = 0, where £_(k) denotes the band with the 
lowest energy. The approximate, normalized ground-state eigenvector has the form 

|*k-) = UVfl_\0) = 4EE Q KM (n t )cl [ -j=\0), (62a) 

where Z k = iV _1//2 exp(— v^/2) and the coefficients are given by 

C iKM (m) = K - J r= r K~ expitk • Rj). (62b) 

Thus, our first task is to evaluate the matrix elements (60b). A somewhat lengthy, but 
straightforward calculation gives the following results 

z 

M£ + = — t k ^2 cos k a + h£lvk(vk — 2g) — 2jgMlv^ 1 exp(— v%.) sinhi£, (63a) 

a=x 

M+~ = -t k (e 2in/3 cos k x + e~ 2in/3 cos k y + cos k z ), (63b) 
M k - = M k ++, M k - i.\l k )*, (63c) 

where j — for the quasi-Holstein model, j = 1/2 for the E®e JT polaron, and 

t k = texp (-*;£) . (64) 

Having obtained the matrix M k , we readily find its eigenvalues £ K (k), 

£ K (k) = -t k E K (k) + hD,v k (v k - 2g) - 2jghQv k 1 exp(-vl) sinhi£, (65) 

where the quantities 



E K (k) = e (k)-K^e 2 (k)+e 2 (k) (k = ±1) (66) 

denote the two e g bands in the absence of the JT coupling and 

e (k) = cos k x + cos k y + cos k z , (67a) 

ei (k) = - (cos k x + cos k y — 2 cos k z ), (67b) 

£2 (k) = -— (cos/ca; — cosA^). (67c) 
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The condition 5£-(k) / 5v^ = then yields a transcendental equation for v^, which may be 
written as 

l + 2jexp(-^)F k 

Uk = 9 i + ^/nm, ' (68a) 

where = E_(k) and 

F k = exp(-t£) - (2u£) _1 sinh^. (68b) 

To assess the range of validity of our variational approach, we shall now first investigate the 
limiting cases of weak and strong coupling, 
(i) weak coupling: g <C 1 

In this coupling range, Eqs. (68) possess the solution 

^ = rP^ + 0(g\ (69a) 

where the adiabaticity parameter 

7 = t/hQ. (69b) 

If this is inserted into the expression for £-(k), Eq. (65), the result is 

£_(k)/^ = - 7 £ k -s 2 -i±^, (70) 

where terms of 0(j 2 ) have been excluded because of the expansion (51). To compare our 
formula with exact analytical results for the Holstein polaron, we set j = and restrict 
ourselves to one dimension. At the T point, Eq. (70) then reduces to 

£ lD (0)/nn = - 1 -- [ ^-. (71a) 

This may now be compared with the result of weak-coupling perturbation theory, 27 which 
is valid for all 7: 

£ _">(0)/Ml = - 7 - 7 =^=. (71b) 
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While these two expressions agree for 7 <C 1, they start to diverge for larger 7, and our 
formula (71a) gradually ceases to be a reasonable upper bound. 26 For the pure JT case 
( 7 = 0, j = 1/2), Eq. (70) reduces to 

£_(k)/Ml = -2g 2 , (72) 

which agrees with the result of perturbation theory. 21 
(ii) strong coupling: g >> 1 

In this case, Eqs. (68) are solved by the expression 

v* = g-j- g + 0(g~ 2 ), (73) 
which, after substitution into Eq. (65), leads to 

S4k)/hn = -g 2 -j. (74a) 
For the Holstein polaron, strong-coupling perturbation theory 27 yields the result 



7 



2 



£.(}L)/Ml=-g 2 -\+) , (74b) 

which agrees with our formula (74a) (for j = 0) in the nonadiabatic limit 7 C 1. In the 
pure JT case (j = 1/2), Eq. (74a) agrees with the strong-coupling expression in Ref. 21. 

This concludes our discussion of the weak and strong coupling limits. Summarizing we 
may state that our variational treatment seems to work reasonably well in the nonadiabatic 
regime 7 < 1, but becomes less reliable for larger 7. 26 The nonadiabatic regime might be 
relevant to the manganites. For, in the doping region considered in this work, the ti g core 
spins form an antiferromagnetic (G-type) spin background, leading to a strong suppression 
of hopping because of the double-exchange mechanism. 9 

Important characteristics of a polaron are its dispersion relation and effective mass. 
Subsequently these properties will be examined for both the quasi-Holstein model (j = 0) 
and the E%e JT case (j = 1/2) by means of a numerical evaluation of Eqs. (68). The 
result of such a calculation for g — 1.5 and 7 = 0.5 is shown in Fig. 1, where the polaron 
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dispersion relations £_(k) for j — and j = 1/2 are depicted along a closed path of the 
cubic Brillouin zone (BZ). For comparison, the lower tight binding e g band (g = 0) is also 
shown in Fig. 1. Although the polaron bands are shifted to lower energies and have a smaller 
width in comparison with the e g band, as expected, the shapes of all three bands are very 
similar. In particular, all extrema of the dispersion curves are located at the same positions. 
Other prominent features of the bands are the extended flat minima between T and X and 
the absolute maxima at the R point. 

Figures 1 and 2 
The polaron effective mass m* is defined by the relation 

J_ J_ / <9 2 £_(k) \ 
m* ~ ¥ \ dk 2 J k=0 ' 

where k is the wave-vector component along some symmetry line of the BZ. Using Eq. (65) 
we obtain the simple result 

m /m* = exp(-vl =0 ), (75) 

where 

J_ = _l ( d2E A 

m h 2 \ dk 2 / k=0 

denotes the inverse effective mass of the lower tight binding e g band. We mention that in 
our derivation of Eq. (75) the k dependence of has been properly taken into account, the 
simplicity of the result being due to a cancellation of all terms involving d 2 vy_/dk 2 . 

In Fig. 2 the mass ratio m /m* is plotted as a function of the coupling strength g. 
Although both polaron masses behave similarly, there is an unexpected crossover where, for 
increasing g, the Holstein polaron starts to acquire a somewhat larger effective mass than 
the JT polaron. The effect is, however, not as dramatic as claimed by Takada. 29 For, in the 
strong-coupling limit, the polaron effective mass ratio mj TP /m* HP tends to the finite limit 
exp(— j) PS 0.607 for j = 1/2, as follows from Eqs. (73) and (75). The results presented 
in Fig. 2 are in remarkably good agreement with recent quantum Monte Carlo data, 14 
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thus confirming our expectation that the proposed variational approach, restricted to the 
nonadiabatic regime, gives a fair account of polaronic properties over the whole coupling 
range. 

VII. SUMMARY 

In this work a detailed account has been given of the analytic properties of the E (g) e JT 
polaron, consisting of a mobile e g electron linearly coupled to the local e g normal vibrations 
of a periodic array of octahedral complexes. The linear JT coupling implies the existence 
of two operators, the angular momentum J and the parity /C, which commute with the JT 
part and are responsible for the twofold degeneracy of all JT eigenvalues. This degeneracy 
is lifted by the anisotropic hopping term, which does not commute with J and /C. The 
most interesting feature of our study is, however, the appearance of a close relationship 
between the JT problem and the Holstein model. Although such a connection has already 
been suspected to exist in the simpler E(g)b JT polaron, 11 it has never been explicitly 
demonstrated. This connection only emerges in a particular representation of the original 
problem, in which the Hamiltonian acquires an explicit dependence on the half-integral 
angular momentum quantum number j and quite naturally decomposes into a Holstein 
term and a residual JT interaction. While the ground state of the JT polaron belongs to 
the sector j = 1/2, the Holstein polaron is formally obtained for the unphysical value j = 0. 
This is the optimal form of the Hamiltonian, which can be achieved by purely analytic 
means, allowing the JT and the Holstein polaron to be treated in a unified framework. 

The Hamiltonian is then subjected to a variational treatment, yielding approximate 
ground-state energies and eigenvectors for both types of polarons. Although the ground- 
state eigenvector is explicitly given by Eqs. (62), its application to the calculation of physical 
properties is relegated to future work. Here we have restricted ourselves to the polaron 
dispersion relations and effective masses. As expected, the polaron bands are shifted to 
lower energies and have a smaller width in comparison with the bare e g band (see Fig. 1), 
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but the shapes of all three bands are very similar. The dependence of the effective masses on 
the coupling strength g is also similar for both polarons (see Fig. 2). There is, however, an 
unexpected crossover where the Holstein polaron starts to acquire a somewhat bigger mass 
than the JT polaron with increasing g. These results are in remarkably good agreement with 
recent quantum Monte Carlo data. 14 This seems to indicate that our variational approach, 
restricted to the nonadiabatic regime for formal reasons, is fairly accurate. The nonadiabatic 
regime might be relevant to the manganites, for in the doping region considered in this work 
the core spins form an antiferromagnetic (G-type) spin background, leading to a strong 
suppression of hopping because of the double-exchange mechanism. 9 

A more realistic model would have to include, at least, the intersite coupling of the 
normal vibrations. This coupling is expected to contribute to the splitting of the degenerate 
JT ground state and to turn the Einstein phonons of the present work into optical phonons. 
Whether these coupling terms will give rise to additional and unexpected effects, remains 
to be seen and will be investigated in future work. 
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APPENDIX A: DIRECT PROOF OF EQUIVALENCE 



In Sec. V the equivalence of the operators A and Aj has been proven in terms of their 
algebraic properties. Here we intend to put forward a more direct proof by showing that 
the operators possess identical matrix elements, provided the states involved are related to 
each other by the mapping prescription (39). We start by defining the basis vectors 

\ik] m + 2n«, m) = d} K \m + 2rij, m) |0/, 0/), (Ala) 
\m] m + 2rii + 1, m + 1) = dJ K \m + 2ui + 1, m + 1) J| |0j, 0/), (Alb) 

where the notation is the same as in Eqs. (28) and k, is fixed by Eq. (29c). The vectors 
(Ala) and (Alb) are then eigenstates of the angular momentum operator J to the same 
eigenvalue j = m + 1/2 and belong to the subspaces Uj~ and U~, respectively, according to 
our definitions introduced at the end of Sec. IV. Given these basis states, the nonvanishing 
matrix elements of the operator A, Eq. (30a), are readily evaluated by means of the easily 
proven relations 

aix\m + 2rii, m) = \Jm + rii \m + 2rii — 1, m — 1), (A2a) 
diz\fn + 2rii, m) = ^Jnl \m + 2n« — 1, m + 1), (A2b) 

H K \m + 2rii, m) = 0, Il_ K |m + 2nj, m) = \m + 2rii, m). (A2c) 

With the help of Eqs. (A2) we then find the expressions 

A\ik; m + 2n«, m) = \/2ni m + 2n» — 1, m + 1), 
A\ik] m + 2rii + 1, m + 1) = \J 2j + 2n^ + 1 \m\ m + 2nj, m), 

and, hence, the only nonvanishing matrix elements of the operator A read: 

(in; m + 2rij — 1, m + l|.A|i/c; m + 2nj, m) = y/2rii, (A3a) 
(in; m + 2n i? m\A\iK] m + 2rii + 1, m + 1) = 2j + 2nj + 1 . (A3b) 
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Our assertion is that the matrix elements of the operator Aj, defined by Eqs. (40), are 
the same as those on the right side of Eqs. (A3), provided the states are chosen as prescribed 
by Eq. (39). To prove our claim, we first introduce the new basis vectors 

\ iK;2ni ) = cl^d=\0), (A4a) 



2n,)! 

\ iK . 2n t + 1) = 4 );> 10), (A4b) 

y/(2ni + l)\ 

where we have adopted the notation of Eqs. (37). Here, in contrast to the vectors (Al), the 
quantum number k needs no longer to be fixed, but may be arbitrarily set equal to 1 or —1. 
Since, by definition, the vectors (A4a) and (A4b) are elements of the subspaces and V~, 
respectively, the mapping prescription (39) requires the following one-to-one correspondence 
to exist between the vectors (Al) and (A4): 

\ik; m + 2rii, m) <-> \ik; 2rii), (A5a) 
\ik; m + 2rii + 1, m + 1) <-> \ik; 2rii + 1), (A5b) 
\m; m + 2rii — 1, m + 1) <-> \ik; 2rii — 1). (A5c) 

By means of these relations the matrix elements (A3) are then mapped onto the following 
expressions 

(in; 2rii — \\Aj\in\ 2rii) = \j2n^ (A6a) 
(i«; 2ni\Aj\iK\ 2n^ + 1) = yj2j + 2rn + 1, (A6b) 

whose validity is explicitly verified with the help of Eqs. (A4) and (40). This proves the 
equivalence of the operators A and Aj, as far as the basis states are concerned. Owing to the 
completeness of these states, however, relations (A3) and (A6) suffice to extend the proof to 
arbitrary vectors of the spaces Uj and Vj, provided these vectors are related to each other 
by Eq. (39). This completes our direct proof of the equivalence of the operators A and Aj. 
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Figure Captions 

Fig. 1. Polaron dispersion curves £_(k)/Ml, Eq. (65), for j — (Holstein polaron) and 
j = 1/2 (E<g)e JT polaron). The lower tight binding e g band (g = 0) is also shown for 
comparison. Symmetry points of the BZ are designated as in Ref. 28. 

Fig. 2. Inverse effective masses of the Holstein polaron (dashed line) and the _E®e JT 
polaron (solid line) as functions of the coupling strenght g. 
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